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We study how to incorporate Mott physics in the BCS-type superconductor, motivated from the 
fact that high T c superconductivity results from a Mott insulator via hole doping. The U(l) slave- 
rotor representation was proposed to take local density fluctuations into account non-perturbatively, 
describing the Mott-Hubbard transition at half filling. Since this decomposition cannot control local 
pairing fluctuations, the U(l) slave-rotor representation does not give a satisfactory treatment for 
charge fluctuations. Extending the U(l) slave-rotor representation, we introduce an SU(2) slave- 
rotor representation to allow not only local density fluctuations but also local pairing excitations. 
We find an SU(2) slave-rotor gauge theory of the Hubbard model in terms of two kinds of collective 
boson excitations associated with density and pairing fluctuations that interact with gapless fermion 
excitations via SU(2) gauge fluctuations. An interesting observation in this effective description is 
that phase fluctuations of fermion pairs arise as SU(2) gauge fluctuations. Thus, fermion-pairing 
excitations can be controlled by dynamics of collective bosons in the SU(2) slave-rotor gauge theory. 
Performing the standard saddle-point analysis based on the SU(2) slave-rotor action, we find an 
interesting phase described by partial freezing of charge fluctuations near the Mott-Hubbard critical 
point, where local density-fluctuation modes are condensed but local pair-excitation modes are 
gapped. Partial freezing of charge fluctuations causes fermion pairing to be incoherent as a result 
of reconciliation of superconductivity and Mott physics. The nature of this non-superconducting 
phase is identified with an anomalous metal due to the presence of incoherent pairing. 

PACS numbers: 71.10.Hf, 71.30.+h, 71.10.-w, 71.10.Fd 



I. INTRODUCTION 

Nature of the doped Mott insulator is one of the central 
interests in modern condensed matter physics, associated 
with the mechanism of high T c superconductivity. Con- 
sidering the high T c phase diagram, an antiferromagnetic 
order in the parent Mott insulator vanishes rapidly via 
hole doping, and a paramagnetic non-superconducting 
phase appears. This non-magnetic state evolves into the 
high T c superconductor, doping holes further. [1] The cen- 
tral issue is the nature of the intermediate non-magnetic 
phase between the antiferromagnetic Mott insulator and 
high T c superconductor. It is important to notice that 
physics of the superconducting state is BCS-like although 
the high T c superconductivity results from the doped 
Mott insulator. [1] In this respect the intermediate phase 
will be determined by competition of Mott physics and 
BCS superconductivity. In this paper we discuss how to 
incorporate the Mott physics into the BCS superconduc- 
tivity. 

To understand the Mott physics in a concrete manner, 
we consider the BCS-Hubbard Hamiltonian 

H = "* E ( C U* + H.c.) 

- £ lAtfle"*"* (4^ n - 4^) - H.c. + i £ | Ay | 2 

<y) (ij) 

+ u ^2 n ^ n H^ (!) 

i 

where Ay is the pairing potential with its amplitude 



| Ay | and phase ay. This Hamiltonian introduces the 
competing nature arising from the density-phase uncer- 
tainty; the J term causes local pairing of electrons while 
the Hubbard-[/ term suppresses local charge fluctua- 
tions, thus breaking phase coherence of electron pairs. 
When local charge fluctuations are strong in the case 
of U < U c with the critical strength U c for the Mott 
transition, phase fluctuations of electron pairs would 
be suppressed, i.e., (e -mij ) = 1 owing to the density- 
phase uncertainty relation. In this case superconductiv- 
ity is expected to appear, described by the BCS model 

L SC = J2ia4<r( d r ~ M) c «x - tJ2(ij)A c L c j^ + H - c -) ~ 
|Ay Kcjjct^— cj^ctj)— if.c. On the other hand, when 
local density fluctuations are suppressed in the case of 
U > U c , the density-phase uncertainty causes (e~ iaij ) = 
0. The resulting non-superconducting phase would be de- 
scribed by the effective Lagrangian L e f j = ^2 icr c\ a (d T — 
- 1 E (y> (cU> +H.c.) - £ (i .> \A i:j \e- ia v (c] T 4 - 

c]; c j-f) — ^' c ' — j? 52n cos (<5 x a )i where the gauge ac- 
tion is introduced to impose (e~ laij ) — 0. This state 
can be regarded as a non-Fermi liquid metal with phase- 
incoherent pairs. [2, 3] 

The problem is how to control local charge fluctu- 
ations. Recently, U(l) slave-rotor representation was 
proposed in order to take local charge fluctuations into 
account non-perturbatively. [4] Actually, the U(l) slave- 
rotor gauge theory of the Hubbard Hamiltonian found 
the bandwidth-control Mott transition from spin liquid to 
Fermi liquid at half filling, ignoring antiferromagnetism 
and superconductivity. [4-6] Using the U(l) slave-rotor 
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decomposition Ci a = e 
expression from Eq. (1) 



{ fieri we obtain the following 



Z = J Df ier D0 i D\A ij \Da ij DL i Dtp i e-tf drL , 



+7Ei^i 



(2) 



where the Hubbard-?/ term is represented as the charge 
and spin channels f/^i^T^a = f Ei^t + n ii - I) 2 - 

f £>it " ^l) 2 + f £>*T + nil) - T Ei 1. and the 
spin channel is not taken into account in this paper. It 
is easy to show that Eq. (2) is exactly the same as Eq. 
(1) with the charge channel for the Hubbard- U term, in- 
tegrating out the ifi and Li fields with = e l6i Ci a . 
In this expression the electron Hilbert space \ci a > is 
reconstructed as the direct product of the boson and 
fermion Hilbert spaces \Li >0|/ i(J > according to the 



decomposition c,, 



-iO, 



fia, where Li represents an 



electron density at site i. It is clear that any decomposi- 
tion method enlarges the original electron Hilbert space, 
thus an appropriate constraint associated with the de- 
composition should be imposed. The Lagrange multi- 
plier field (fi expresses the U(l) slave- rotor constraint 
Li = ^2, a fia — 1, implying that the fermion and boson 
Hilbert spaces are not independent, thus the two opera- 
tors fi a and e l0i also. Then, e~ l9i is identified with an 
annihilation operator of an electron charge owing to the 
constraint Li = f i(7 fia — 1 an d the canonical relation 
[Li,9j] = —iSij imposed by —iLid T 9i. 

Integrating out the density field Li and performing 
the Hubbard-Stratonovich (HS) transformation to de- 
compose the "correlated" hopping term in the following 

way 



/ dr{-tJ2(fl^e-^f ia +H.c.)\ 

Jo {ij)° 

= J DctijD(3ijexp - J dr|t^(a*j/3jj + H.c.) 

° {ij) 

-t E (fl°tjfj° + H - c -) f E^/v-** + h.c.)}] , 



{ij)v 



we obtain the U(l) slave-rotor gauge Lagrangian with 



d — wave pairing 

Z = J Df itT D6 i Da ij DtpiDci j e-tf d - rL , 

A 2 

L = L + Lf + Lg, L = ^(ta/?+— ), 

{ij) 

L f = E fl^r - *I>E flf* !] 

ia i <j 

-taY,{fle- iCii fj°+H.c.) 

{ij)° 

-a ^ ftj . e -**« (44 - 44) - h.c, 



{ij) 

Le =rJ E^ ~ Vl)2 ~ U zZ^" ieinj<: i0j + Hx ^ 

(3) 



{ij) 



where the hopping parameters are represented as a y - = 
ae lCij and = f3e lCij with amplitudes a, (i and phase 
, and d — wave pairing-symmetry is assumed in qj of 
| Ay | = Aq-ij with the gauge transformation a y - — » a y - + 
6*i + 6>j . The unidentified hopping and pairing amplitudes 
can be determined self-consistently in the saddle-point 
analysis, a = \( e -^e ie i +H.c.)\, = |(£ ff flf ja +H.c.)\ 

and A^ = |(44-44)|. 

This effective Lagrangian is quite appealing. Start- 
ing from the BCS-Hubbard effective model, we ex- 
tract dynamics of collective charge fluctuations from the 
Hubbard-{7 term. As a result, the fermion sector corre- 
sponds to the BCS description with a renormalized band- 
width Da while dynamics of collective density excitations 
is described by the boson-Hubbard-type model with an 
effective chemical potential in the saddle-point approx- 
imation (cjj = 0). Remember that the boson-Hubbard 
model is the prototype for the quantum phase transition 
associated with charge fluctuations, describing the gen- 
uine Mott transition without symmetry breaking via con- 
densation. The quantum transition indeed occurs when 
D(3/U ~ 1 with the half bandwidth D, where D(3 is an 
effective bandwidth for the boson field. One can show 
that the hopping parameter (3 decreases as U increases, 
resulting in the bandwidth-control Mott-Hubbard tran- 
sition. Actually, the previous slave-rotor studies showed 
this transition without superconductivity (A = 0) at half 
filling. [4-6] 

However, this treatment does not take into account 
phase fluctuations e _lay of fermion pairs carefully al- 
though local density fluctuations are governed by U(l) 
rotor excitations e l6i . In other words, the U(l) slave- 
rotor representation cannot give any condition for pair- 
ing excitations. Considering the mathematical structure 
of Eq. (3), one can find two kinds of gauge excitations, 
c^ and a-ij associated with local density fluctuations and 
pairing excitations. Density-gauge excitations (cjj) can 
be controlled by U(l) rotor excitations e \ When local 
density fluctuations are suppressed in the large U limit, 
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boson excitations are gapped (Mott insulator), causing 
gapless density-gauge fluctuations. On the other hand, 
when local density fluctuations become strong in the 
small U limit, the U(l) rotor variables get condensed, 
making density-gauge fluctuations massive due to the 
Andcrson-Higgs mechanism. However, there are no such 
boson excitations corresponding to the U(l) density-rotor 
variable for pair-gauge fluctuations (a^ ) . 

This discussion motivates us to introduce new boson 
excitations associated with local pairing fluctuations, al- 
lowing us to control phase fluctuations of pairing exci- 
tations (pair-gauge fluctuations). Assume the presence 
of such boson excitations. We can estimate that both 
density and pairing fluctuations would be gapped due to 
the density-phase uncertainty in the large U limit. On 
the other hand, both collective excitations become con- 
densed in the small U limit. What happens near the Mott 
critical point? It should be noted that local pair excita- 
tions are different from local density fluctuations. Thus, 
there is no reason for both excitations to be condensed 
at the same time. Can there exist an intermediate phase, 
where only one kind of boson excitations is condensed. If 
possible, what is the nature of this intermediate phase? 

In this paper we extend the U(l) slave-rotor formula- 
tion, allowing local pairing fluctuations. Then, the col- 
lective boson field is expressed as an SU(2) matrix field 

{ Zii -zl \ 

Uiaa* — 1 + ] j involved with both density (z^) 
\ M z l\ J 

and pairing (zi\) fluctuations. [7] This leads us to con- 
struct an SU(2) slave-rotor gauge theory of the Hubbard 
model in terms of the collective density- and pairing- 
fluctuation modes interacting with gapless fermion exci- 
tations via SU(2) gauge fluctuations. [8] Compared with 
the U(l) slave-rotor gauge theory [Eq. (3)], z^ corre- 
sponds to e l9i roughly speaking while Zi± is newly intro- 
duced to control pairing-gauge fluctuations (phase fluc- 
tuations of electron pairs). Since the SU(2) slave-rotor 
gauge theory admits two kinds of collective bosons, an 
intermediate phase can be naturally allowed between the 
spin liquid Mott insulator {(z^) — 0, (zn) — 0) and 
Fermi liquid metal ((z^) ^ 0, (z^) ^ 0), characterized 
by softening of the density- fluctuation modes ((z^) ^ 0) 
as the Fermi liquid, but gapping of the pair-excitation 
modes ((zjj) = 0) as the spin liquid. Remember that the 
U(l) slave-rotor gauge theory contains only the density- 
fluctuation modes (z^), thus allowing the two phases 
only, the spin liquid Mott insulator ((z^) = 0) and Fermi 
liquid metal ((zjf) ^ 0). 

Performing a mean-field analysis based on an SU(2) 
slave-rotor effective action, we find the intermediate 
phase indeed away from half filling, described by par- 
tial freezing of charge fluctuations ((z^) ^ 0, (zi±) = 0) 
near the Mott-Hubbard critical point. We reveal that 
this intermediate phase originates from emergence of 
a pseudospin-dependent chemical potential due to hole 
doping that causes SU(2) pseudospin symmetry breaking 
to discriminate the local pair-excitation modes (zjjj from 



the local density-fluctuation modes (z^). Condensation 
of density-fluctuation modes makes density-gauge excita- 
tions gapped due to the Anderson-Higgs mechanism. On 
the other hand, pairing-gauge fluctuations remain gap- 
less because pairing-excitation modes are gapped. As a 
result, we find an effective U(l) gauge theory for this 
intermediate phase, where phase fluctuations of fermion 
pairs are described by pairing-gauge fluctuations. We 
discuss physics of this non-superconducting phase, iden- 
tified with a non-Fermi liquid metal due to the presence 
of phase-incoherent pairs. 



II. SU(2) SLAVE-ROTOR THEORY OF THE 
HUBBARD MODEL 

A. SU(2) slave-rotor representation in the path 
integral formulation 

We start from the Hubbard Hamiltonian 

H = -* E + H.C.) + y]T ct T C iTC t iCa ,(4) 

(ij)a i 

where u/t is a parameter determining various phases of 
this model with hole concentration. The local interaction 
term can be decomposed into the pairing and density 
channels in the following way 

3u t t u Y^ t t 

yL^t^i = 2 E C *T C W C *T 

i i 

+lE(E^W-i) 2 + |(E4^-i). 

1(7 (T 

Performing the HS transformation for the pairing- and 
density-interaction channels, we find an effective La- 
grangian in the Nambu-spinor representation 

Z=[ D[^,*f, ^^ile-Zo"^, 

l = E 4(Pri - - tj2(4^j + h.c.) 

1 (ij) 

i 

+^-X^ 2+ ^ 2+ ^- (5) 



Here ipi is the Nambu spinor, given by ipi = ( t ) • ®i 

V c i.i ) 

and 3>f are the real and imaginary parts of the on-site 
pairing-order parameter respectively, and tpi is an effec- 
tive density-potential, ji is an electron chemical poten- 
tial, renormalized by local interactions as fi = Hb — u/2, 
where /j,f, is the bare chemical potential. Introducing a 
pseudospin vector = (<f>f , $j,(fi), one can express Eq. 
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(5) in a compact form 



(8) reads 



Z = J D[1> i ,ft i ]e-ti dTL , 

» (ij) 

-i^Vj^-r^ + ^-^tr^.f) 2 . (6) 



Integrating over the pseudospin field fij , Eq. (6) recovers 
the Hubbard model Eq. (4). Note that the U(l) slave- 
rotor representation allows only the density-interaction 
channel including the t 3 matrix. 

As discussed in the introduction, we disintegrate bare 
electrons into collective excitations and rcnormalizcd 
electrons in the following way 



ipi = e 



BliTl — I<p2i1"2— «<p3iT 3 



where the two component spinor rji 



(7) 



can be 



considered to express renormalizcd electrons, and the 
SU(2) matrix field Ui = cxp[z J2k=i ^ki^k] collective 
bosons. Here exp[i^ij"n] (exp[?0 2 iT2]) can be interpreted 
as a creation operator of an electron pair since it mixes 
a particle with a hole. If one takes only the (f>u phase 



field in the Ui matrix field, i.e., Ui 



m , he finds 



r] i+ = cos (piCtf+i sin (fiicl^ and r/|_ = i sin fac^ +cos0jcjj 
from r\i = e z ^ uTl tpi. Considering the <p2i phase field, one 
obtains = cos^jCif+sin^jct, and r]J_ — — sin^jCjj + 

cos <pic\, in the same way. On the other hand, exp[i<j) 3 iT 3 ] 
is identified with a creation operator of an electron 
charge, corresponding to the rotor variable in the U(l) 
slave-rotor representation. 

Inserting the SU(2) slave-rotor decomposition Eq. (7) 
into the effective Lagrangian Eq. (6), we obtain 



Z = J D[rH,Ui,{U]e-ti*TL, 

L = Y, vKdrl + UAUj - nUiTsUj)^ 

i 

-t^ivjUinUjvj+H.c.) 

(ij) 

-1 vluSi ■ f)Uj Vi + i- tr(^ • ff. (8) 



Performing the unitary transformation fij • f — > U\ (Oj • 
T)Ui that makes the integration- measure invariant, Eq. 



L = J2 4(drl + Uid T U} - nUinUj)^ 

i 

-tJ^ivlUmUfrj + H.c.) 

(ij) 

-i Vt@i ■ r)Vi + ^E tr ^' f ) 2 - ( 9 ) 



Shifting f2j • f as Oj • f — iUid T U\ + i\iUiT 3 U\ that also 
makes the integration- measure invariant, we find the 
SU(2) slave-rotor representation of the Hubbard model 
as an extended version of its U(l) slave- rotor representa- 
tion 

Z = J D[rH,Ui,(U]e-^ dTL , 

L = J2vl(d T I - i(U ■ r) Vi -t^ivlUiTsUjrjj + H.c.) 

i (ij) 
+ ^ £ tr (~ iU i d rU} + (U ■ T + inUiT 3 Ujf. (10) 



B. Relation between the path integral formulation 
and canonical quantization 

One cautious person may suspect the above deriva- 
tion because any decomposition method should always 
impose its associated constraint, but the path integral 
derivation does not seem to include an appropriate con- 
straint. However, this guess is not correct because the 
above derivation imposes its constraint indeed. [6, 9] 

We start from the easy-axis anisotropy with Oj = 
(0,0, (fii) and Ui = e z ^ 3T3 in order to check whether Eq. 
(10) recovers the U(l) slave-rotor Lagrangian [Eq. (3) 
with A = 0] in this limit. Inserting this easy- axis repre- 
sentation into Eq. (10), and performing the HS transfor- 
mation for the last time-fluctuation term, we obtain 

Z = J D[r]i,(f) 3 i,ipi,Li]exp - J dT^Y^4 d rVi 



, T 3 e- l4 '^ T3 T lj + H.c.) 



(ij) 



£Q L ? ~ iL i( d T<f>3i - if*) + i<Pi(Li - Vi^m))} 



(11) 

where r\i and e~ 4< ^ 3iT3 carry spin and charge degrees of 
freedom, respectively. In this expression the electron 
Hilbert space \tpi > is decomposed into the direct product 
of the boson and fcrmion Hilbert spaces \Li \ Vi > ac ~ 
cording to the decomposition ipi — e~ t ^' iiT3 r]i 1 where Li 
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represents an electron density at site i. Since this de- 
composition enlarges the original electron Hilbcrt space, 
the U(l) slave-rotor constraint Li = r]jT 3 r]i is introduced 
to reduce the enlarged Hilbert space into the original 
electron one. Then, e~^ 3iT3 is identified with an an- 
nihilation operator of an electron charge in the Nambu- 
spinor representation owing to the constraint = rjjrsrji 
and the canonical relation [Li,(f> 3 j] — —i5ij imposed by 
—iLid T 4>3i- Actually, Eq. (11) is the starting point in 
the canonical quantization of the U(l) slave-rotor repre- 
sentation for the Hubbard model, Eq. (3) with A = as 
discussed in the introduction. [4, 6, 9] In this respect the 
above expression connects the canonical approach with 
the path integral representation. 

In the same way as the above, performing the HS trans- 
formation for the last time-fluctuation term in Eq. (10) 



we find an effective Lagrangian of the Hubbard model 

Z = J D[ m ,U l ,n i ,E lJ ,F l0 ]e~^ dTL , 
L = Lq + L n + Ljj, 
L =tJ2tr(F ij Ej j +E ij F} j ), 



L v = E vlidrl - i&i ■ r)Vi -t^ivjFijrjj + H.c.) 



-t^2tT(U}E} j U i T 3 + H.c), 

(ij) 



(13) 



Am 



where Eij and Fij arc HS matrix fields associated with 
hopping of r\i fermions and Ui bosons, respectively. 
We make an ansatz for the hopping matrix fields as 



^trjuMf + iMi(-iUid T Uj + ■ f + i(j,UiT 3 U})}, E ij ~ £cxp[m y - • t\t 3 , F l3 w Fexp[m y - • t]t 3 (14) 



we obtain the following expression for the SU(2) slave- 
rotor representation 

L = J2 4drVi -tJ2(vlUir 3 U] Vj + H.c.) 



(ij) 



+ 



£ [tr{«M? + iMi(-iUid T Uj + ifiUiTsll})} 



• (tr[Mif| - 4^Vi) 



(12) 



where E and F are longitudinal modes (amplitudes) of 
the hopping parameters, and a,j their transverse modes 
(phase fluctuations) , considered to be spatial components 
of SU(2) gauge fields. The reason why we introduce the 
t 3 matrix is that the fermion sector L v should recover the 
original electron Lagrangian [Eq. (6)] as the slave-rotor 
representation does. 

Inserting Eq. (14) into Eq. (13), we reach the SU(2) 
slave-rotor gauge theory of the Hubbard model for the 
metal-insulator transition 



where Mi is an SU(2) pseudospin matrix. In this 
SU(2) case the electron Hilbert space is represented as 
\ipi >— |Mi >0|?7i > according to the decomposition 
ipi = Ujrji. The Lagrange multiplier field fij confirms 
the local constraint tr [Mif] = "qlrrji. Using Mj = N ■ f 

where Ni corresponds to the pseudospin-density-wave or- 
der parameter, one can show that this is nothing but the 
standard relation Ni — \ri\fr]i as a natural extension 
of the U(l) slave-rotor constraint. Lfj is identified with 
an annihilation operator of an electron pseudospin in the 
same way as the U(l) case. 



C. SU(2) slave-rotor gauge theory of the Hubbard 
model 



Using the HS transformation for the correlated hop- 
ping term in Eq. (10) 



-t{r}\ a Uiai3T 3l 3 1 U} l5 'n jS + H.c.) 
— ► t Fij a $E^ Sa + Eij a $F^ Sa 
-(ViaFijaSVjS + U ial3 T Wl U] lS El iSoi ) - H.c. 



Z = J D[rH,Ui,(U,3ij\e-^ dTL , 
L = L v + Lu + 4tY,EF, 



(ij) 



1 (ij) 

L v = ^J2tr(-iUid T Uj + & ■ t + inUmU}) 2 

i 

-tEY,^(U]T 3 e-^- ? U t T 3 + H.c), 

(ij) 

where the unknown parameters such as the amplitudes 
of the hopping parameters E, F and the SU(2) pseu- 
dospin order parameter Oj should be determined by the 
following self-consistent mean-field equations 

4£ = (rjl^rij + H.c), 4F = {trOj}T 3 U t T 3 + H.c.)), 

= iuirilrrH) - ^(tr[f(-iUid T Uj +i»UiT 3 U})}}(16) 

as will be discussed later. 

Compared with the U(l) slave-rotor effective action 
Eq. (3), the SU(2) slave-rotor effective action Eq. (15) 
exhibits more fruitful physics. Dynamics of collective 



(15) 
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charge fluctuations can be extracted from interacting 
electrons in the strong coupling approach as the U(l) 
case, but the structure of their dynamics is much richer. 
The effective boson action is represented as the nonlin- 
ear a model- type (its CP 1 representation) with an effec- 
tive chemical potential (the time component of an SU(2) 
gauge field) in the saddle-point approximation ignoring 
SU(2) gauge fluctuations while that in the U(l) slave- 
rotor representation is governed by the boson-Hubbard 
model-type. The presence of additional collective charge 
fluctuations opens the possibility of new phases, as will 
be discussed in this paper. 

Considering the fermion sector of the SU(2) slave-rotor 
representation, one cautious person can find the possibil- 
ity of superconductivity in the Hubbard model. See the 
kinetic energy term of the fermion sector in Eq. (15). 
Because the SU(2) gauge matrix Wij = exp[— ictij ■ r] 
has nonzero off-diagonal components, fermion pairing as 
the seed for superconductivity is naturally allowed. Lo- 
cal repulsive interactions (the Hubbard-C/ term) appear 
to be local density and pairing fluctuations of the SU(2) 
matrix field Ui, governed by the nonlinear a model- type. 
These charge fluctuations generate SU(2) gauge fluctu- 
ations and couple to them, one of which corresponds to 
phase fluctuations of fermion pairs, observed in the ki- 
netic energy term of the fermion sector. Since the SU(2) 
gauge fluctuations are controlled via the collective charge 
fluctuations, fermion-pairing fluctuations can be man- 
aged by the collective boson-dynamics, varying the cou- 
pling constant u/t and hole concentration 8. In summary, 
the SU(2) slave-rotor representation not only reveals the 
possibility of superconductivity in the Hubbard model, 
but also controls phase fluctuations of fermion pairs. 



D. U(l) charge-rotor gauge theory of the Hubbard 
model 

Since the structure of the SU(2) slave-rotor action [Eq. 
(15)] is complex to analyze, it is necessary to consider its 
easy-axis limit f2j-r = ifiTs, Ui = exp[z</>3jT3], and dij-f = 
a 3y T 3- Then, Eq. (15) is reduced to the effective U(l) 
gauge Lagrangian in the slave-rotor representation[4-6, 
8, 9] 

L V = E Viidrl - i<PiT3)Vi 
i 

-tF^nle^^Wj + H.c), 

(ij) 

L <t> = 2~ E^ 3 * - (fit- V) 2 

i 

-2fcEy^cos(<fej - <j> 3i - a 3i j). (17) 

Based on this slave-rotor effective Lagrangian, Florens 
and Georges performed a saddle-point analysis at half fill- 
ing, and found a coherent-incoherent transition of the foi 



fields, identifying this transition with the Mott-Hubbard 
transition from a spin liquid Mott insulator to a Fermi 
liquid metal. [4] For our later discussion we perform the 
mean- field analysis of the U(l) slave-rotor theory. 

For the saddle-point analysis one can resort to the large 
N generalization for the boson sector replacing e l<t>3i with 

Lx = h J2^ X l -b- Vi]Xl)(drX ia + [ M - Vi ]X i<7 ) 



(18) 



where A^ is a Lagrange multiplier field to impose the 
uni- modular constraint for the rotor field X i(7 , and ipi 
is substituted for iipi. Based on the mean- field ansatz of 
<Pi = ipo, «Aj = A, and a 3i j = 0, we find the mean-field 
free energy 



(iv + fi- ip ) 2 



+ \ + Ee*} 

Y,{DEF - fxS - A), (19) 



13^^ 1 2u 



where ui (v) is the Matsubara frequency for fermions 

(bosons) with inverse temperature (3 = 1/T, and e^ X ^ = 
— 2i(cos fc x +cos k y ) is the fermion (boson) dispersion that 
originates from the electron dispersion. 8 is hole concen- 
tration. 

Minimizing the free energy with respect to E, F, <p , A, 
and /z, we find the self-consistent saddle-point equations 
for the mean-field parameters 

DE=[ D deD { e)W ( , 
J-D M 2 -(^o-Fe)2' 

™~ /><«>! s 
/ 



M +\ + Ee 



a nn 1 V" ~ Fe 



u 



-8= \ deD{e)-Y^ 

J-D U ^- fc 



1 



iv+jx—i£o 



±t-^y+x + Et 



■ (20) 



Here J2k ^ s re pl acc d with J D D deD(e), where D is the 
half bandwidth and D(e) is the density of states for an 
electron band. 

Performing the Matsubara frequency summation and 
energy integration with the constant density of states 
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D(e) = l/(2D) in Eq. (20), we obtain the algebraic equa- 
tions for the mean-field parameters E, F, ip a , A, and /i 
at zero temperature 



1 



( <Po V 

\df) 



(2A - DE)y/2u{\ + DE) 



3(DE) 2 
-(2A - Ee_)^j2u(X + Ee_ 

1 (M-^o) 2 
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1 



A 



i = 



DE 2u 
v/2u(A + DE) - y / 2u{X + Ee-) 



A 



DE + DE 



where e_ is given by e_ = 



DE 
1 (M-1o) 2 



2m 



(21) 
From 

__ DFS and 

/i = -DFS - ^2u[X - DE(l - 8)}. Inserting these ex- 
pressions into the above equations, we obtain 



(p-yo) 

2u 



- A 



the third and fifth equations we find ipo 



E = 
F 



1-S 2 



1 



'i{DEf 



(2A - DE)^2u(X + DE) 



-[2A + DE(1 - 8)} v / 2u[X- DE(l-S)} , 
_ v /2u(A + DE) - x /2u[X - DE(l - S)} 



DE 



(22) 



Condensation of the <j)^i bosons occurs when their excita- 
tion gap closes, given by /J. c — f0c = or A c — DE c (l — S) = 
that determines the Mott-Hubbard critical point 



Uc _ 1 / 1 - S 2 \ 
D ~ 4K2-5 ) 



(23) 



with F c = (1 — 2<5)/3 in the mean-field approxima- 
tion. This means that the spin liquid Mott insulator 
((e^ 3i ) — 0) appears in u/D > u c /D, and the Fermi 
liquid metal ((e** 3 *) ^ 0) arises in u/D < u c /D. It is 
important to notice that since the slave-rotor decompo- 
sition is meaningful in the case of E > and F > 0, the 
critical value of F c = (1 — 25) /3 indicates that this rep- 
resentation is valid when hole concentration is relatively 
small, here S < 1/2. 



E. Spin liquid Mott insulator away from half filling 

One would be surprised at this result that the slave- 
rotor theory allows the spin liquid Mott insulator away 
from half filling, in contrast with the common belief that 
hole doping to a Mott insulator is expected to cause a 
metal. Hole doping to a Mott insulator shifts the chemi- 
cal potential from the middle of the charge gap generated 
by the energy cost (u) for double occupancy to the top of 



the lower Hubbard band. [10] Since the density of states 
at the top of the lower Hubbard band is nonzero, the 
resulting state is expected to be a metallic phase. Actu- 
ally, the t — J model studies find a metallic phase at all 
fillings, even arbitrarily close to half filling. [11] 

Recently, Choy and Phillips claimed that this common 
belief may not be right in the Hubbard model. Doped 
Mott insulators can be insulators. [10] They demonstrated 
that hole localization can obtain because the chemical po- 
tential lies in a pseudogap which has vanishing density of 
states at zero temperature. The pseudogap in the doped 
Mott insulator results from short-range antifcrromag- 
netic correlations corresponding to the nearest-neighbor 
singlet-triplet splitting. They showed that the pseudogap 
vanishes without the triplet contribution which lies above 
the chemical potential, claiming that the presence of the 
upper Hubbard band is crucial for the pseudogap. More 
fundamentally, they proposed that physics is sensitive to 
the order of limits of u — > oo and L — > oo, where L is the 
system size. They suggested the relevance length scale 
£do for the pseudogap, where ^ Q represents the average 
distance between doubly occupied sites. They claimed 
that the order of u — > oo and L — > oo results in £<j > L, 
metallic transport while ^ Q < L and localization obtains 
in the reverse order of limits, provided that nh^ a < L, 
rih = x(L/a) 2 the number of holes with lattice spacing a. 
Furthermore, considering the scaling form Z ~ [ J -( t / u ) T 
with p > for the one-hole quasiparticle weight Z, they 
addressed that the discrepancy between the t — J and 
Hubbard results implies the non-commutativity of the or- 
der of limits. In the t — J model (no double occupancy) 
corresponding to u — > oo and L — > oo, Z remains finite. 
But, Z vanishes in the reverse order of limits (Hubbard 
model). 

The slave-rotor theory admits double occupancy. In 
addition, it was demonstrated that the exchange energy 
scale J ~ t 2 /u is indeed generated at the level of one- 
particle properties, where it cuts the divergence of the 
effective mass at the Mott transition driven by u. [4] Actu- 
ally, the single-particle spectral weight is enhanced near 
the chemical potential with vanishing density of states 
at the chemical potential, decreasing dimensionality from 
infinite dimensions near the Mott critical point at half fill- 
ing, which is associated with the exchange energy scale. [4] 

To see the mechanism for the spin liquid Mott insula- 
tor away from half filling in the slave-rotor theory more 
concretely we consider the effective chemical potential 
in the boson sector of Eq. (18) carefully. It is impor- 
tant to notice that the effective chemical potential for 
the rotor condensation is given by \i — tp instead of 
cpo or a. The point is that the effectiv e chemical po- 
tential [i — ipo = —y/2u[X — DE(1 — S)] away from half 
filling shows essentiall y the same b ehavior as the chem- 
ical potential fi = —y/2u[X — DE] with (fo — at half 
filling although fi = -DFS - y/2u[X - DE(l - 5)] and 
fj£>o = —DFS away from half filling have doping depen- 
dencies proportional to hole concentration. [12] Consider- 
ing physics of the slave-rotor variable, this behavior of 
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the effective chemical potential can be understood. The 
U(l) slave-rotor variable raises or lowers the local density 
of fermions. Thus, hole concentration has nothing to do 
with the density of slave-rotor bosons. This is in contrast 
with the slave-boson representation of the t — J model, 
where hole concentration is exactly the same as the den- 
sity of slave-bosons in the saddle-point approximation. 
In the t — J model the single occupancy constraint of 
fivfio- + b\bi = 1 with spinon f ia and holon 6, results in 
ifiafia) = 1 — S and (b\bi) = S, thus causing the con- 
densation of holons due to the holon chemical potential 
as soon as holes are doped. On the other hand, there 
is no such constraint in the slave-rotor representation, 
thus the density of slave-rotor bosons does not follow 
hole concentration. Hole doping changes the density of 
fermions only ((rjlT^rii) — —5), reflected in the doping de- 
pendence of the effective chemical potential <po = —DF5 
for fermions. The effect of hole doping on the boson 
dynamics is just to modify the bandwidth from DE to 
DE(1 — $). As a result, the slave-rotor excitations can 
be gapped away from half filling. 

Compared with the work of Choy and Phillips, [10] the 
slave-rotor theory seems to have some similarities. The 
effective "charge" chemical potential fi — ipo, associated 
with softening of local charge fluctuations, can lie in the 
gap away from half filling in the large u/t limit, similar to 
the half- filled case, while the effective "spin" chemical po- 
tential (fo shows the conventional behavior as the Fermi 
liquid or the fcrmion dynamics in the U(l) slave-boson 
theory of the t — J model. Considering that the slave- 
rotor description can capture local antiferromagnetic cor- 
relations of J, allowing double occupancy seems to play a 
crucial role for the emergence of the Mott insulator away 
from half filling. 



F. U(l) pair-rotor gauge theory of the Hubbard 
model 

Now we consider the easy-plane limit Oj • r = ^fri, 
Ui = exp[z</>ijTi], and • f = aujT\. Then, we find 
another effective U(l) gauge Lagrangian[13] 

i 

-tF^vl^'^r^+H.c.), 

(ij) 

i 

-2tE^2cos{(t>ij+(t>u-auj)- (24) 

(ij) 

We note that Eq. (24) can be reduced to Eq. (6) 
with Oj • t = <&f n after the gauge transformation of 
$f -> $f + d T 4>u, anj -> a U j + 4>u + 4>ij, and Eq. 
(7) without <p2i and <fei are utilized. If the T3 matrix is 



not used in Eq. (14), the hopping term in vanishes, 
and Eq. (6) cannot be recovered from Eq. (24). Ignor- 
ing gauge fluctuations, and replacing $f with $0 as the 
mean-field approximation, we obtain the on-site pairing 
order parameter given by $0 = iu(r)jT\rii) + (d T <pu)- It 
turns out to be zero because double occupancy costs too 
much energy. Thus, there is no phase transition in the 
fcrmion Lagrangian at half filling as the case of the easy- 
axis limit. 

To examine the boson Lagrangian in the mean-field 
level, we resort to the large N generalization as the 
U(l) slave- rotor representation. Introducing the N- 
componcnt rotor field Y^, we rewrite in Eq. (24) 
as 

Ly = Yu E(^ r -)(^ r -) - tE E + Y i° Y *) 
+*]>>(]ri^i 2 -i), (25) 

i a 

where is a Lagrange multiplier field for the uni- 
modular constraint. If we represent Eq. (25) in terms of 
Yia = Ria + ilia (Y^ = R ia - we obtain the follow- 
ing expression for the mean-field Lagrangian of collective 
pair excitations 

LY = ^Y t [{d T R ia f + {d T I ia f\ 

-2tE J2 {Ri*Ri* + iM + * E x & R i« + J -) - 

(ij)cr i & 



One can find that Eq. (26) is exactly the same as Eq. 
(18) if we rewrite Eq. (18) in terms of = Ri„ + ili a . 
Thus, the mean-field analysis in the previous slave-rotor 
theory can be directly applied to Eq. (25). This leads 
us to conclude that both the (j)^ and <j>u fields are simul- 
taneously incoherent in (u/t) > (u/t) , and coherent in 
(u/t) < (u/t)o at half filling, where (u/t)o is the critical 
value for the Mott transition obtained in Eq. (23). This 
originates from ip = at half filling and u c = at the 
Mott-Hubbard critical point, making Eq. (18) exactly 
the same as in Eq. (24). Fundamentally, the rea- 
son why both fields should be coherent simultaneously is 
the presence of the SU(2) pseudospin symmetry at half 
filling. The slave-rotor action should be symmetric (in- 
variant) under the transformation <j>u < — ► </> 3i at half 
filling in the mean-field approximation. 

III. MEAN-FIELD ANALYSIS OF THE SU(2) 
SLAVE-ROTOR THEORY 

So far, we discussed that only two phases are expected 
to appear at half filling in the context of the SU(2) slave- 
rotor theory, corresponding to the spin liquid Mott insu- 
lator with (e^) = 0, (e**) = and the Fermi liquid 
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metal with (e^ 3i ) ^ 0, (e**) ^ 0, respectively. There is 
no intermediate phase because the presence of the SU(2) 
pseudospin symmetry makes both phase fields simulta- 
neously coherent or incoherent. Away from half filling 
an intermediate phase, characterized by condensation of 
only one kind of bosons, is expected to arise because the 
SU(2) pseudospin symmetry is broken explicitly due to 
hole doping, reflected in the chemical potential term. 

For the saddle-point analysis we consider the following 
mean-field Lagrangian from Eq. (15), where SU(2) gauge 
fluctuations are ignored (e Wij T — > I), and the easy-axis 
anisotropy for the SU(2) pseudospin order parameter is 
naturally chosen (Oj • f — > — iipor 3 ), 

L = L v + L v + 4i^£F, 

(ij) 

1 (ij) 

L u = ]-J2 tr(-iUid r U} - vpon + ijJiUiTzU}) 2 



4m 

-tE^riUjrsUm + H.c). 

(ij) 



(27) 



The first assumption to ignore the SU(2) gauge fluctu- 
ations is the simplest mean-field ansatz. The second 
assumption of the easy-axis anisotropy for the SU(2) 
pseudospin order parameter can be justified from the 
self-consistent mean-field condition = iu^lrrn) — 
\{tr[f{-iUid T U} + ifj,UiT 3 U})]). The first term in this 
expression becomes zero for n and t 2 because double 
occupancy costs too much energy. The second term 
also vanishes for t\ and r 2 because this term contains 
pseudospin-flipping terms in the representation Ui — 



-z 1 

H 

t 



of the SU(2) slave-rotor matrix, but the 



boson Lagrangian L\j has no pseudospin-flipping terms 
as will be seen below. Particulary, the second mean-field 
ansatz has important physical implication for the SU(2) 
slave-rotor theory that the SU(2) pseudospin symmetry 
is broken away from half filling. 

Representing Eq. (27) in terms of the z ia field, we 
obtain the following expression 



L = L n + Lu + 2t^2EF, 

(ij) 

L n = Yl Viid-r 1 - ¥>or 3 )??i - tFj^inl^ + H - c -), 

i (ij) 
^ J2( drZ i* - & ~ Ol P0\zt){drZi 

(azlz ja + H.c.) + 



Lagrangian Eq. (28) seems to be similar to the U(l) 
charge-rotor mean-field Lagrangian Eq. (17) with Eq. 
(18), there is a crucial difference between them. The 
most important ingredient in Eq. (28) is the pscudospin- 
dependent chemical potential crip in the boson sector. 
We show that this results in the condensation of only 
one kind of bosons, allowing an intermediate phase away 
from half filling. Since the renormalized dispersion of 
the Zi a bosons depends on the pseudospin component, 
the lowest energy of the z^ bosons lies at momentum 
(0, 0) while that of the zn bosons appears at momentum 
(tt,tt). There are no pseudospin-flipping terms in the bo- 
son Lagrangian, thus the pseudospin-flipping terms such 
as (zj^d T Zii) and (zj^Zii) appearing in the self-consistent 

condition for the SU(2) pseudospin order parameter fij 
should be zero, justifying our assumption of the easy-axis 
anisotropy for the order parameter. 



Integrating out the rji and z i{ 
mean-field free energy 



fields, we obtain the 



F„ = -5^h[(»)'-( w -ftB ! l 



+ J2(DEF - fiS - A). 



(iu + fi- aifo) 2 



+ A + Eae%] 



(29) 



Since the Zii bosons have their energy minima at mo- 
mentum Q = (it, 7r), wc shift the momentum k of the zn 
boson to k + Q. Then, the mean-field free energy is given 
by 



FMF = -- p Y J H(^f-(^-Felf] 



k.uj 



{iv + /i - cT(p ) 2 



k,v cr 

+ YSPEF-n&-\), 



+ A + Eel 



(30) 



-tE 



2u 
E 



where the a symbol disappears in the boson dispersion. 
Compared with the U(l) charge-rotor free energy Eq. 
(19), Eq. (30) allows the pseudospin-dependent chemical 
potential aipo in the boson free energy. The charge chem- 
ical potentials appear to be /j, — ip for the z^ field and 
[/j, — aipo]zia) ^^~fo f° r the field with the spin chemical potential tp a 
in the SU(2) slave-rotor theory while the charge chemi- 
cal potential is given by \i — f n for the e*^ 3i field with 



1), (28) 



where IE and t 3 are replaced with E and a — ±, re- 
spectively. Although the SU(2) slave-rotor mean-field 



the spin chemical potential ipo in the U(l) charge-rotor 
theory. 

Minimizing the free energy with respect to E, F, ipo, A, 
and /i, we find the self-consistent saddle-point equations 
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for all mean-field parameters 



H at zero temperature 



J-D u ^H 2 -(^-f f ) 2 ' 

DF = — f D deD(e)^-yy — % , 

J-D P V a ~ (w+ »- u aipo) + X + Ee 

J-D { 'eW-^o-ft) 2 



1 



A + E'e 



±ra^ + A + Ee 



(31) 



One can see from the denominator of the boson propa- 
gator that condensation of the z ia bosons occurs when 
their excitation gap closes, given by — [i c + ip 0c = for 
the Zi-[ bosons at momentum (0, 0) and — \i c — ip 0c = 
for the bosons at momentum (jr, 7r), respectively. For 
both bosons to be condensed simultaneously, pL c = and 
ipoc = should be satisfied. Actually, this happens at half 
filling where both bosons have exactly the same transi- 
tion point given by A c = E C D as a result of the SU(2) 
pseudospin symmetry Away from half filling the effec- 
tive chemical potential (po becomes nonzero, causing the 
SU(2) pseudospin symmetry breaking not spontaneously 
but explicitly because it is pseudospin-dependent. As 
a result, the bosons should be condensed before the 
condensation of the Zii bosons occurs. Only one kind 
of bosons can be condensed in the intermediate range of 
u/t, corresponding to the intermediate phase mentioned 
before. 



Performing the Matsubara frequency summation and 
energy integration in Eq. (31), we find the following ex- 
pressions for the mean-field parameters E, F, ipg, A, and 



E = 



1 



i-(-5*-V 

\DFJ 



6(DEy 



(2A - DE)y/2u{\ + DE) 



-(2A - £e_) v /2w(A + Et_) 
1 



-(2A - Ee + )y/2u(\ + Ee+) 

^o(^ + puj = 0, 

y/2u(X + DE) - v /2ti(A + Ee_) 



1 = 



2DE 



y/2u(X + DE) - v /2u(X + Ee+) 



+ 
5 = 1 



A 



2DE 
1 



DE 2DE 



(M- Vo) 2 (m + ^o) 2 



2u 

(ft±(Po) 2 



2u 



(32) 



2u 



- A 



One can 



where e± is given by e± = 

see the validity of these equations performing the limit 
of 6 — > 0. We obtain cp = at half filling from the 
third equation. This causes [i = — y / 2u(X — DE) in the 
fifth equation. Using the expressions for e± with /j, = 
— y/2u(X — DE) and ip a = 0, one can find e± = —D, 
reproducing the mean-field equations for the half-filled 
case. [12] 

As soon as holes are doped, th e chemical potential /i 
jumps from (i = —y/2u(X — DE) to fj, = —(E/F)u. De- 
creasing the Hubbard interactions with this chemical po- 
tential, one can find the Mott-Hubbard transition point, 
where the z^ bosons begin to be condensed. The z^ con- 
densation occurs when A c + E c e_ = is satisfied. From 
the analytic expressions in Eq. (32) with e±, we obtain 
the following conditions 



E, 



1 



fAd 



u c E c r A c _ u c fE c \ 2 ' 

d f c Y d d\fJ y 



(33) 



determining the critical value u c /D for the Mott- 
Hubbard transition in the SU(2) slave-rotor theory away 
from half filling. We note that the above equations for 
the Mott critical point do not recover the half-filled case 
owing to the chemical potential jump. From these equa- 
tions we find u c /D w 0.170 at 5 = 0.010, u c /D w 0.171 at 
5 = 0.050, and u c /D w 0.173 at 5 = 0.100. The value of 
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FL (-t) 18 NFl(-t) 28 SL 



<zt>^0 <zt>^0 <zt>=0 u/t 
<Z|>=jfcO <z ; >=0 <z ; >=0 

FIG. 1: A phase diagram in the SU(2) slave-rotor theory of 
the Hubbard model away from half filling 



u c j D for the Mott critical point increases as hole concen- 
tration becomes larger, consistent with our expectation. 
Remember that this Mott critical point is defined by the 
condensation of the z^ bosons while the bosons are 
gapped owing to — u c — Lp 0c > 0. We identify this Mott 
critical point with (u/t) 2 s- 

There is another critical point associated with the con- 
densation of the Zii bosons with (z%-\) ^ 0. It is clear 
that both the z^ and z^ bosons get condensed when 
A c + E c e- = (-fi c + tpoc = 0) and A c + E c e + = 
(— /i c — ipoc = 0) are satisfied, i.e., /j, c = and <poc = 0. 
Then, we obtain A c = (1 — 5)DE C from the last equation 
in Eq. (32). Inserting this expression with E c = 1/2 
into the fourth equation in Eq. (32), we find the second 
quantum critical point 



1 



D 4,(2-6)' 



(34) 



The actual values are u c /D w 0.127 at S = 0.010, u c /D w 

0. 128 at 5 = 0.050, and u c /D w 0.132 at 6 = 0.100. 
At half filling Eq. (34) recovers the Mott critical point 
u c /D = 0.125 of the U(l) slave-rotor theory exactly. We 
define this critical value as (u/t) is. Note that (u/t) is < 
(u/t) 2 s is satisfied as expected. 

In summary, there is an intermediate region (u/t)is < 
u/t < (u/t) 2 s away from half filling where only the z^ 
bosons are condensed, i.e., (z^) ^ and (z^) = 0. The 
phase in u/t < (u/t)is is characterized by (z^) ^ 
and (zii) ^ 0, identified with the Fermi liquid metal 
while the phase in u/t > (u/t) 2 s is given by (z^) = 
and (zii) — 0, named as the spin liquid Mott insulator. 
We emphasize again that the presence of the interme- 
diate phase in (u/t)is < u/t < (u/t) 2 s results from the 
pseudospin-dependent chemical potential aipo due to hole 
doping, causing the pseudospin SU(2) symmetry break- 
ing. The mean-field phase diagram is summarized in Fig. 

1, where "FL", "NFL", and "SL" represent Fermi liquid, 
non-Fermi liquid, and spin liquid, respectively. Notice 
that the intermediate phase is identified with the non- 
Fermi liquid metal in the phase diagram Fig. 1. In the 
following we investigate the nature of this intermediate 
phase named as the non-Fermi liquid metal. 

One cautious person may suspect that the intermediate 
phase is an artifact of gauge degrees of freedom. In other 
words, performing gauge rotations to this ground state 
results in the same phase as the Fermi liquid state. How- 



ever, this guess is partially correct only. It is important 
to discriminate the pseudospin SU(2) symmetry from the 
SU(2) gauge symmetry. As discussed before, hole doping 
causes the easy-axis anisotropy, breaking the SU(2) pseu- 
dospin symmetry explicitly. However, the SU(2) gauge 
symmetry cannot be broken. [14] Unfortunately, the inter- 
mediate phase breaks the SU(2) gauge symmetry. Is this 
phase a gauge artifact? Not necessarily. This situation is 
exactly the same as the characterization of superconduc- 
tivity. Usually, we characterize superconductivity as the 
condensed phase of Higgs bosons (Cooper pairs). This 
identification necessarily breaks the U(l) gauge symme- 
try, thus cannot be right in a rigorous manner. In this 
respect there is no local order parameter to break the 
U(l) gauge symmetry. However, we should admit that 
this characterization is useful. The way to understand 
this identification is to break the gauge symmetry ex- 
plicitly, fixing one gauge. Since the gauge symmetry is 
explicitly broken via gauge fixing, one can define the lo- 
cal order parameter such as Cooper pair bosons. Now we 
can understand the intermediate phase in the same way. 
Consider gauge fixing first. Then, measure the conden- 
sation of SU(2) rotor bosons. This is a physically appeal- 
ing way to discriminate various phases. In the following 
section we show that the intermediate state is indeed dif- 
ferent from the Fermi liquid phase owing to the presence 
of incoherent fermion pairs. 



IV. NATURE OF THE INTERMEDIATE PHASE 



To reveal the nature of the intermediate phase, it is 
necessary to find an effective Lagrangian. Low energy 
boson excitations are considered to be 



-yT — xe 
sjx 



(35) 



for this ground state. Here x = \(z^)\ 2 is the conden- 
sation amplitude, not determined self-consistently in the 
previous saddle-point analysis owing to complexity. 0i 
is a dynamic variable to guarantee partial freezing of 
charge fluctuations, i.e., (zi{) = ^Jl — x(e l6i ) = 0. The 
ansatz of Eq. (35) satisfies the uni-modular constraint 
5^ CT \zi a \ 2 — 1. Since the objective of this section is to 
find a physical picture for this intermediate phase, we do 
not try to determine the condensation amplitude x, and 
just assume its presence. 

The above ansatz for the SU(2) matrix field leads the 
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time-fluctuation term of the boson sector to be 

i 

^ ^[(1 - x)(d T 6i + fl zl - itpo - in) 2 



hopping term is written as 
-tEj^triUfaWijUiVi) 

ij 

= -tEY\x{X^ + 4 T ) - (1 - X){e-^X}^ 



—2y/x{l — x)(Cl x i cos Qi + flyi sin 9i)(d T 9i — 2i/j.) 

+2(n 2 xi + n 2 j + 2x(n zl - i<po + ifif 



+e^X in e- i9i } - ^x{l-x){x\ n e^ +X loi e- 10 



(36) 



where represents the time component of the SU(2) 
gauge field around its zero mean-value. The 2x(Q zi — 
i(fo + iii) 2 term allows us to replace fl Z i with iip — ifi in 
the low energy limit. This corresponds to the Anderson- 
Higgs mechanism, as will be also seen in the boson hop- 
ping term. Inserting Q zi = iip — i[i into Eq. (36), we 
obtain the following expression for the right-hand-side 
(RHS) in Eq. (36) 



RHS 
l-x 



+X m e-»*+X} jl e»>} 



(39) 



2u 



^ (d T 0i - 2ijj, - J [n xi cos 9i + n yi sin 9^ 



i 

One can see Q x j and £l y i fluctuations gapped from 

i 

= I E (J 1 - x cos2 9 A + 1 1 ~ x sin2 e ^li 

i 

\ [Q xi cos 9i - n yi sin6y 2 J, 



The first term in the RHS makes the Ay-j bosons con- 
densed. This is nothing but the Andcrson-Higgs mecha- 
nism because the condensation of the zy bosons makes 
the gauge fluctuations in the Ay-p fields gapped. This 
corresponds to the fact that the <\>zi condensation causes 
gapped a^ij excitations in the U(l) slave-rotor gauge the- 
ory. Considering the uni-modular constraint, we can set 
A i: ,j = y/y and Xy-| = y/1 — ye laij as Zy — y/x and 
Zii = y/l — xe ; of Eq. (35), where y is the condensation 
amplitude. Then, the boson hopping term is obtained to 
be 

-tE^^iUjrsWijUm) 

ij 

= -tEYjpty/V - 2(1 - aOv^cos^j - 9j) 



u 



-2y/x{l - x)(l - y){cos(9i - ay) + cos(9 3 - ay )} 
- -tEY\^x^y-2{l-x) s /ycos{9 i -9 j ) 

ij 

—2yjx{l — x){\ — y){cos(9i — 9j — ay) + cos ay} 

(40) 

where we shift the gauge field the 
last line. 

An important issue arises due to the cos ay potential 
term. One may claim that this term allows us to ignore 
the gauge fluctuations at low energies because it causes 
ay excitations gapped. However, this problem does not 
seem to be so simple owing to the stiffness parameter in 
the cos term. When the condensation amplitudes x and 
y of the Zi-\ anc ^ X ij1 bosons respectively are close to 1, 
the stiffness parameter of the cos potential becomes very 
small. This may cause gauge fluctuations gaplcss. One 
can argue that the cos potential would be always rele- 
vant in two space and one time dimensions [(2 + 1)D] 
in the renormalization group sense. However, this claim 
can be applied to the conventional sine-Gordon model 
in (2 + 1)D.[15] The present problem is more complex 
since gauge fluctuations also couple to fermion excita- 
tions, resulting in screening and dissipation in the gauge 
dynamics. 

First, we discuss the case when the cos ay term is 
relevant. Then, we shift ay as 5y + n due to the 
minus sign of the cos potential, and set ay + tt = 
with the constraint | Ay 1 1 2 + 1 Ay J 2 = 1. Then, the boson ay — 9j + it = owing to the relevance of the cos 



thus safely ignored in the low energy limit. As a result, 
we find the low energy part for the time-fluctuation term 



l-x 
2u 



Y,{d T 9 l -2in) 2 . 



(37) 



To obtain the low energy sector for the boson hopping 
term, we use the following representation for the SU(2) 
gauge matrix Wy 



X H\ ~ X ijl 



(38) 
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term. This is an important constraint for gauge fluc- 
tuations. Because we consider the intermediate mean- 
field phase characterized by (z;j) = \/l ~~ x(e l6i ) = 
and (z^) = y/x, gauge fluctuations cause (e laij ) = 
that results from the constraint (e WiJ ) = (e*w -7r )). As 
a result, the SU(2) gauge matrix Wij is found to be 



- ye 



with (e lo «) = 0. 



Note that this SU(2) gauge matrix is consistent with the 
mean-field analysis in Eq. (28) since this expression is 
reduced to the unit matrix in the mean-field level owing 
to (e taij ) = 0. We emphasize that the off-diagonal com- 
ponents of the gauge matrix can arise in the intermediate 
phase beyond the mean-field approximation, associated 
with fermion pairing. 

Using this gauge matrix, we find the following expres- 
sion for the fermion hopping term 



(41) 



The key feature in the fermion hopping term is that 
fermion pairing is dynamically generated beyond the 
mean-field approximation due to the off-diagonal com- 
ponents of the SU(2) gauge matrix although this fermion 
pairing vanishes in the mean-field fashion owing to 
(e* ay ) = 0. This fermion hopping term identifies the 
off-diagonal gauge fluctuations with phase fluctuations 
of fermion pairs. Since the rji fermion carries an electric 
charge owing to the condensation of the z^ bosons, we 
interpret these dynamically generated fermion pairs as 
preformed Cooper pairs owing to strong phase fluctua- 
tions of the fermion pairs. An interesting point is that 
the preformed pairs arise from the kinetic energy term 
instead of the potential term in the Hubbard model. 

The low energy effective Lagrangian for this phase is 
found to be 

icr (ij) a 
i 

-2tEj2lVx(±-x)(±-y) - (1 - x)y/y\ cos^ - 6j), 



(ij) 



(42) 



where we use Cl z i = iipo — ip in the time-fluctuation term 
for the rji fcrmions, obtained in the boson hopping term. 
Remember the gauge constraint dij+n = aij — 9j+n = 0. 
One can suspect this derivation because the pairing sym- 
metry is s — wave instead of d — wave. In fact, the 
pairing symmetry depends on the sign of the gauge ma- 
trix Wij. Thus, one can obtain the d — wave fermion 



pairing from the kinetic energy term considering Xij± — 
—y/\ — y<,ije laii , where is + when j = i ± x, and — 
when j = i ± y. In this case the low energy effective 
Lagrangian becomes 

L v = Yl vL( d r - v)Vicr - tFy/y (vLvja + H.c.) 

icr (ij)a 

-tF^T^Y.^'^^A- - - H - c > 

(ij) 



-2tE J2 [ftj Vx(l-x)(l-y) -(l-x)^y\ cos(0j - 

(ij) 

(43) 

It is not easy to determine which pairing symmetry will 
appear without more sophisticated analysis. However, 
one can see that if he compares the boson Lagrangian in 
Eq. (43) with that in Eq. (42), the stiffness parameter 
in Eq. (43) is larger than that in Eq. (42). In the case 
of d — wave pairing we find p x =tE\ y/ x{l ~ x)(l — y) + 
(1 ~x)y/y\ and p y = tE\yJx(\ - x) (1 - y) - (1 - x)^/y\, 
where p x (y) 1S the stiffness parameter in the x (y) di- 
rection. On the other hand, we obtain p x = p y = 
tE\y/x(l — x)(\ — y) — (1 — x)^/y\ in the case of s — wave 
pairing. Because the stiffness parameter in the d — wave 
case is larger than that in the s — wave case, we ex- 
pect that d — wave pairing may be more favorable than 
s — wave pairing. 

Note that if the condensation amplitude x is close to 1, 
the assumption of (z^) = \f\ — x{e l9i ) — is consistent 
with the small stiffness parameter in Eq. (43). In this 
case the resulting fermion Lagrangian becomes with the 
electromagnetic vector potential Aij 

i<7 

-tF^yY^le^Vjc+H.c.) 

(ij)v 

- t Fy/i^Y,^ e ~ iaii ^iA- - - H - c - 



(ij) 



\ Y COS ( d X fl )' 
$ n 



(44) 



where the last gauge action is introduced to impose the 
condition (e~ laij ) — with an internal gauge charge g of 
the rn fermion. 

Eq. (44) is the main result of this section. Based on 
this effective Lagrangian for the intermediate phase, we 
discuss its physical implication. If quantum corrections 
due to gauge fluctuations a^j are ignored as the mean- 
field approximation, the effective fermion Lagrangian is 
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reduced to 

i<7 

-tFrf^inle^mr + Hx.). (45) 

This seems to coincide with the Fermi liquid metal. Al- 
lowing the density-rotor variable only as the U(l) slave- 
rotor representation, the Fermi liquid metal appears 
when the charge-rotor bosons are condensed. [4] Actually, 
the charge-rotor condensation makes the fermion field rji 
couple to the electromagnetic field Aij. However, the 
presence of the pair-rotor variable is expected to alter 
this physical picture since the r]i a fermions are not elec- 
trons but fractionalized ones due to gapped pair-rotor 
excitations. To see the existence of coherent electron ex- 
citations, we consider the electron spectral function given 
by the convolution integral between the r\i„ and z i(T prop- 
agators 

G e m{ij,TT f ) = (T T [c TiT c{ jV ,]) 

= (TAUlvirVl'Ujr'hl) ~ x(T T [m T+ T$ T , + ]) 

+ (1 - x)(T T [e-»*e»*>])(T T [Tjl T _T, jT ,-]). (46) 

There exist coherent electron excitations in the mean- 
field fashion, resulting from the condensation of the 
bosons. In this respect one can say that the interme- 
diate phase corresponds to the Fermi liquid state. But, 
it should be noted that this result obtains in the saddle- 
point approximation ignoring gauge fluctuations Oy . The 
presence of gauge fluctuations alters this picture com- 
pletely. Notice how the Fermi liquid renormalizes as a 
result of strong interactions. When x = 1, this renor- 
malized Fermi liquid metal becomes identical with the 
Fermi gas. Away from x = 1 we find the quasiparticle 
weight Z e i ~ x, reduced by the presence of gapped pair- 
rotor excitations, where the remaining portion 1 — x is 
transferred into incoherent backgrounds. 

Strong gauge fluctuations corresponding to phase fluc- 
tuations of fermion pairs should be allowed. Then, wc 
expect that this phase can be identified with a non- 
Fermi liquid metal owing to pairing fluctuations. The 
presence of pairing-phase fluctuations is the hall mark to 
discriminate this phase from the Fermi liquid metal, rep- 
resented as U(l) pair-gauge excitations. It is well known 
that the presence of long-range gauge interactions can 
result in non- Fermi liquid physics. [16] When gauge fluc- 
tuations are minimally coupled to gapless fermions, thus 
screening of gauge fluctuations occurs via particle-hole 
excitations, the effective U(l) gauge theory is charac- 
terized by the dynamical critical exponent z = 3. In 
this case the imaginary part of the fermion self-energy 
is given by w 2 / 3 at the Fermi surface, implying that its 
real part also has the same frequency dependence via 
the Kramer's Kronig relation, thus giving rise to a non- 
Fermi liquid behavior. [16] The 7 coefficient of the specific 
heat is proportional to —InT in three spatial dimensions, 



and T -1 / 3 in two dimensions. [17] The dc conductivity 
is proportional to in three dimensions and T~ 4 / 3 

in two dimensions. [18] However, this non- Fermi liquid 
physics cannot be applied to the present case because 
screening of gauge fluctuations arises in particle-particle 
excitations instead of the particle-hole channel. Unfortu- 
nately, we don't know the role of pair-gauge fluctuations 
in the non- Fermi liquid physics at present. This should 
be investigated near future. 

On the other hand, if the 9i fluctuations are sup- 
pressed, i.e., (zii) = \/l — x(e l6i ) 7^ 0, the corresponding 
gauge fluctuations would be also suppressed owing to the 
gauge constraint. In this case the resulting effective La- 
grangian is obtained to be 

i<r 

-tFy/vY,{nle iA *>r, ja +H.c.) 

-tF^T^Y.^Ui- - VU}+) - H.c (47) 

This is nothing but the effective Lagrangian for the 
d—wave BCS superconductivity. Thus, the d—wave BCS 
superconductivity can result from softening of pairing ex- 
citations in the SU(2) slave-rotor gauge theory. The tran- 
sition nature from the anomalous metal with incoherent 
pairing to the d — wave superconductor is beyond the 
scope of this paper. 

The SU(2) slave-rotor gauge theory can find the Fermi 
liquid, superconductivity non-Fermi liquid, and spin liq- 
uid in principle. The Fermi liquid is described by < 
(zii) = \fx < 1 with any (zn) and (X^) = y/y = 1, 
Xiji — (Wij = I). The first condition results in the 
coherent electron-quasiparticle weight given by Eq. (46). 
The second condition indicates the absence of both di- 
agonal and off-diagonal fluctuations of the SU(2) gauge 
matrix fields, where the absence of diagonal gauge fluc- 
tuations originates from the Anderson-Higgs mechanism 
due to the Zi\ condensation while the complete suppres- 
sion of off-diagonal gauge fluctuations should be deter- 
mined by the self-consistent mean-field analysis, not per- 
formed in this paper. In our mean-field analysis Eq. (28) 
both the diagonal and off-diagonal bosons are condensed, 
i.e. (z ia ) 7^ 0, and the complete suppression of the off- 
diagonal components in the SU(2) gauge matrix is as- 
sumed, i.e., Wij = I. 

Superconductivity appears when < (z^) = \fx < 1, 
< (zti) = y/l ~x~ < 1 and < (Xiji) = VV < 1. 
< (Xiji) = VI ~ V — 1- The presence of the off- 
diagonal Zii condensation is necessary for the off-diagonal 
gauge bosons (X^i) to be condensed, causing coherent 
Cooper pairs for superconductivity given by Eq. (47). 
Unfortunately, we could not see superconductivity in our 
mean-field analysis because the off-diagonal components 
in the SU(2) gauge matrix are not allowed from the start. 
It is necessary to perform more sophisticated mean-field 
analysis allowing the off-diagonal gauge components. 
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The anomalous metal arises when < (z^) = \fx < 1, 
(z tl ) = y/ T=x( e**) = and < (Xy T ) = y/y < 1, 
= ^/l - y(e lai i) = 0. One important thing is that 
we cannot see the non-Fermi liquid phase in the mean- 
field approximation ignoring gauge fluctuations. To ob- 
tain the anomalous metallic phase, off-diagonal compo- 
nents should be allowed in the SU(2) gauge matrix Wy, 
meaning that the amplitudes of the off-diagonal compo- 
nents should be nonzero, and the expectation values of 
the off-diagonal components should be zero due to their 
strong gauge fluctuations, not captured in the present 
saddle-point analysis. What we have shown in our mean- 
field analysis is that the intermediate phase characterized 
by partial freezing of charge fluctuations, (z^) ^ and 
(zii) = 0, can appear near the Mott-Hubbard critical 
point. We see that this intermediate phase exhibits the 
Fermi liquid physics in the saddle-point approximation. 
However, the pair-gauge fluctuations (the off-diagonal 
gauge components) should be allowed, modifying the 
Fermi liquid physics completely, (z^) — \/l — x(e lSi ) = 
causes (^yj.) = — y(e laij ) = self-consistently. As 
a result, this intermediate phase is described by Eq. (44) 
with incoherent pairing fluctuations. 

The spin liquid Mott insulator is given by (z^) = 0, 
(zii) = and (X^) = 0, (X^i) = 0. This is consistent 
with the mean-field analysis Eq. (28) except -Xy-j = 0, 
implying that this phase is identified with the U(l) spin 
liquid Mott insulator in the context of our mean-field 
analysis. On the other hand, if one admits fluctuations 
of the SU(2) gauge matrix fields, this insulating phase is 
interpreted as the SU(2) spin liquid Mott insulator due to 
SU(2) gauge fluctuations. In the case of U(l) gauge fluc- 
tuations it was claimed that the U(l) spin liquid Mott in- 
sulator can be stable against confinement resulting from 
instanton excitations when there exist sufficiently large 
flavors of gapless fermion excitations. [19, 20] However, 
in the case of SU(2) gauge fluctuations the SU(2) spin 
liquid Mott insulator would be unstable against con- 
finement due to the SU(2) gauge fluctuations [21, 22] al- 
though there are no reliable calculations for this confine- 
ment problem owing to its complexity If the confine- 
ment is realized at low energies, the SU(2) spin liquid 
is expected to turn into an antiferromagnetic insulator 
due to particle- hole confinement. An interesting possi- 
bility can arise that deconfined fermion excitations usu- 
ally called spinons may appear at high energies owing to 
the asymptotic freedom[21] of the SU(2) gauge theory. 
This may explain the unidentified broad spin spectrum 
in the antiferromagnetic phase of high T c cuprates, ob- 
served at high energies in the inelastic neutron scattering 
measurements. [23] 

So far, we discussed the case when the cos ay term 
is relevant. Now we consider the case when the cos ay 
term is irrelevant. Shifting ay as + 7r due to the mi- 
nus sign of the cos potential, we obtain the SU(2) gauge 

matrix W ti w ( ia .. ^ T_7 ^f ). An im- 

portant difference from the previous case is that there 



is no gauge constraint for ay and 9j because the cos ay- 
term is assumed to be irrelevant. In this and 
6j are all strongly fluctuating. The low energy effective 
Lagrangian is obtained to be 

L v = vhidr - V)Vi* - tFy/y (vhvjv + H.c.) 

itr (ij)cr 

-tFv/I - ^ <,yc {r,\ + r,]_ - r,\_r,] + ) - H.c, 

(ij) 

2u 

i 

-ZtE^Jfiij \J x(l -x)(l- y) cos(9 l - 6j - oy) 

ij 

-(l-x)Vycos(^-^-)]- (48) 

In the intermediate phase {e lSi ) = the fermion La- 
grangian L v in Eq. (48) recovers Eq. (44), thus the 
anomalous metal is allowed even in this case. On the 
other hand, when the zn bosons become condensed, 
the a,ij fluctuations would be suppressed due to the 
Andcrson-Higgs mechanism, thus ay- also, resulting in 
superconductivity described by Eq. (47). However, this 
conclusion is in contrast with the assumption that the 
cosa y term is irrelevant. Thus, to be self-consistent be- 
tween the assumption and conclusion, the cos ay term 
should be relevant. 



V. DISCUSSION AND SUMMARY 

In this paper we try to answer to the question how to 
incorporate Mott physics (charge fluctuations) in BCS- 
type superconductivity. Particulary, the objective is to 
construct an effective theory controlling pairing fluctu- 
ations in the BCS-like effective model. Since the U(l) 
slave-rotor gauge theory allows local density fluctuations 
only, this decomposition cannot control pairing excita- 
tions. Remember that there arc two kinds of gauge 
fluctuations in the U(l) slave-rotor gauge theory with 
fermion pairing [Eq. (3)], one of which corresponds to 
phase fluctuations of fermion pairs. Although the U(l) 
charge-rotor variable governs dynamics of density-gauge 
fluctuations, such boson excitations to control pair-gauge 
fluctuations do not exist in the effective U(l) gauge the- 
ory [Eq. (3)]. In this paper we find the missing col- 
lective charge fluctuations represented as the pair-rotor 
variable in the SU(2) slave-rotor gauge theory. The core 
of the SU(2) slave-rotor gauge theory is that two kinds 
of collective boson excitations are introduced to control 
two kinds of gauge fluctuations. One of the boson exci- 
tations corresponds to the density-rotor variable of the 
U(l) slave-rotor representation, and the other is asso- 
ciated with on-site pairing fluctuations, called the pair- 
rotor variable. One of the gauge fluctuations is inter- 
preted as the density-gauge field of the U(l) slave-rotor 
representation, and the other is identified with phase de- 
grees of freedom for fermion pairs, called the pair-gauge 
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field. Although the SU(2) gauge structure is complex, we 
demonstrate that the density-rotor variable controls the 
density-gauge field while the pair-rotor variable governs 
the pair-gauge field. Since dynamics of collective charge 
fluctuations can be controlled via the local-interaction 
strength u/t, pairing fluctuations are naturally handled. 

We perform the saddle-point analysis based on the 
SU(2) slave-rotor gauge theory. The crucial point in 
this analysis is that the SU(2) pseudospin symmetry is 
explicitly broken via hole doping, causing the easy-axis 
anisotropy. The easy-axis anisotropy of the SU(2) pseu- 
dospin order parameter is reflected as the pseudospin- 
dependent effective chemical potential for the SU(2) 
slave-rotor variable, allowing that only local density fluc- 
tuations are softened while local pairing fluctuations re- 
main gapped in the intermediate parameter range away 
from half filling. Since density-rotor excitations are con- 
densed, this phase exhibits the Fermi liquid physics in the 
saddle-point approximation corresponding to its nonzero 
quasiparticle weight proportional to the condensation 
amplitude of the density-rotor bosons. However, we see 
that this intermediate phase differs from the Fermi liq- 
uid state beyond the mean-field approximation allowing 
gauge fluctuations. Because pair-rotor excitations are 
gapped in this phase, their corresponding pair-gauge ex- 
citations are strongly fluctuating. Remember that the 
pair-gauge fluctuations can be identified with phase fluc- 
tuations of fermion pairs. Thus, incoherent fermion pair- 
ing does exist in this phase. The presence of preformed 
fermion pairs discriminates this phase from the Fermi 
liquid state completely [Eq. (44)]. 

It is valuable to discuss intuitively why local charge 
fluctuations can be partially frozen away from half fill- 
ing. At half filling an on-site density fluctuation should 
induce an on-site pair excitation. This is the origin why 
the density- and pair-rotor bosons should be coherent 
or incoherent simultaneously Away from half filling the 
on-site density fluctuation need not give rise to the on- 
site pair fluctuation because the density fluctuations can 
occur between |0) and |1) sites, but the on-site pair fluc- 
tuations should appear between |1) and |1) sites, where 
|0) is an empty site and |1) is one-electron site. As a re- 
sult, away from half filling the density-rotor bosons can 
be condensed against the on-site Coulomb repulsion us- 
ing the |0) and |1) sites while the pair-rotor bosons are 
gapped due to the on-site Coulomb repulsion in the in- 
termediate parameter range. 

Several important issues remain open. The SU(2) 
slave-rotor gauge theory opens the possibility of super- 
conductivity in the Hubbard model. Introducing the 



SU(2) slave-rotor matrix causes the SU(2) gauge matrix 
naturally. Since the SU(2) gauge matrix allows its off- 
diagonal components, off-diagonal fermion pairing can 
appear naturally in the kinetic energy term. Unfortu- 
nately, we did not perform a mean-field analysis to de- 
termine the amplitudes of the off-diagonal components in 
the SU(2) gauge matrix owing to its complex structure. 
More sophisticated self-consistent saddle-point analysis 
should be performed near future for superconductivity 
in the Hubbard model. Another problem is how to un- 
derstand the spin liquid Mott insulator away from half 
filling. Although the slave-rotor gauge theory admits 
the presence of the Mott insulator away from half fill- 
ing owing to its special structure allowing double occu- 
pancy, more extensive numerical simulations should be 
performed for the Hubbard model. To understand the na- 
ture of the non-Fermi liquid phase with preformed pairs 
such as transport and thermodynamics, it is indispens- 
able to find how to treat pair-gauge fluctuations. 

The stability of the non-Fermi liquid metal with in- 
coherent fermion pairing against disorder is an interest- 
ing open question for this phase to be a genuine metal- 
lic phase in two dimensions. The present author in- 
vestigated the role of disorder in the two dimensional 
fermion system with long-range gauge interactions, where 
gauge fluctuations couple to charge currents instead of 
pair currents as the present case. [6, 24] Long range in- 
teractions are shown to make the fermion system stable 
against weak disorder even in two dimensions because 
the gauge interactions let the fermions lie in a critical 
phase, causing that the critical fermions feel their effec- 
tive dimension higher than two owing to their anoma- 
lous critical exponents. Remember that there exists the 
localization-delocalization transition in the weak disor- 
der limit above two dimensions. Criticality can protect 
fermions from localization due to disorder in the weak 
disorder limit. [6, 24, 25] However, the present case should 
be addressed more thoroughly because gauge fluctuations 
couple to pair currents instead of charge currents. This 
problem may be related with the anomalous metallic be- 
havior in two dimensions. [26, 27] 
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